A van der Waals density functional study of chloroform and bromoform on graphene 
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A calculational study of the trihalomethanes chloroform (CHCI3) and bromoform (CHBra) ad- 
sorbed on graphene is presented. The study uses the van der Waals density functional method 
vdW-DF to obtain adsorption energies and adsorption structures for these molecules of environ- 
mental concern. In this study chloroform is found to adsorb with the H atom pointing away from 
graphene, with adsorption energy 357 meV (34.4 kj/mol). For bromoform the calculated adsorption 
energy is 404 meV (39.0 kj/mol). The corrugation of graphene as seen by chloroform is small, the 
difference in adsorption energy along the graphene plane is less than 6 meV. 

PACS numbers: 31.15.E-,71.15.Mb,71.15.Nc 



I. INTRODUCTION 

Trihalomethane (THM) molecules are small molecules 
that are similar to methane (CH4) but with three of the H 
atoms replaced by halogens (F, CI, Br, I, At). The most 
welknown THM is trichloromethane (CHCI3), also known 
as chloroform. THMs are of environmental concern as 
they are toxic to human health [TJ |2] . The human body 
adsorbs THMs by inhalation and by passage through the 
skin, but the main contribution to human exposure arises 
from the consumption of chlorinated drinking water [3] . 

Of the THMs, chloroform is found in the highest con- 
centration in the environment. Chlorine used for water 
disinfection reacts with organic material in the water, 
forming a number of THMs as byproducts: mainly chlo- 
roform, but also THMs with one or more Br atoms. The 
THMs toxity motivates a search for an effective process 
for selective extraction. Carbon materials, such as acti- 
vated carbon or carbon nanotubes are used or have been 
suggested for the use in adsorbing filters for removing 
THMs from the drinking water after the disinfection, but 
before intake 0] . 

We here study how a chloroform molecule adsorbs on 
the simplest of carbon materials, graphene. By use of 
density functional theory (DFT) calculations we deter- 
mine the energy gained at adsorption and compare with 
the adsorption energies of similar molecules, like methane 
and tribromomethane (CHBrs), also called bromoform. 
For these calculations we apply the first-principles van 
der Waals (vdW) density-functional method vdW-DF 

mm- 

Chloroform on carbon materials has previously been 
studied in a few studies by use of theory. DFT has been 
used for chloroform on benzene in a study employing the 
vdW-DF method [7] and for a study of chloroform on 
carbon nanotubes with use of the local density approx- 
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imation (LDA) 0J. For experiments, there is a century 
long tradition of studies of chloroform because it was fre- 
quently used as a solvent and as an anasthetic. However, 
adsorption studies on carbon materials, provinding ad- 
sorption (or desorption) energies, are more recent [HH1H]. 

The outline of the rest of the paper is as follows: In 
Section II we describe the computational method and 
our system of chloroform and graphene. In Section III 
we present our results and discussions, and Section IV 
contains our summary. 



II. METHOD OF COMPUTATION 

THMs are molecules with a central C atom and four 
other atoms surrounding the C atom approximately 
evenly distributed. Of these, one atom is an H atom 
and the three other atoms are halogens. In this paper we 
analyze the adsorption on graphene of chloroform and 
bromoform, for which the three halogen atoms are CI or 
Br atoms, respectively. 

We use DFT with the vdW-DF method 016] to deter- 
mine the adsorption energy and atomic structure. Our 
calculations are carried out fully self-consistently [5J. 
We use the DFT code GPAW Qlj] with vdW-DF 
in a Fast-Fourier- Transform implementation [11] , The 
GPAW code is an all-electron DFT code based on pro- 
jector augmented waves [T2] (PAW). 

Figure [T] illustrates the adsorbed chloroform molecule 
on graphene, and the periodically repeated orthorhombic 
unit cells used in our calculations. We show the adsorp- 
tion configuration that has the H atom pointing away 
from graphene. A previous vdW-DF study [7] of chloro- 
form with benzene discussed the C-H/7T interaction in a 
dimer with the H atom of chloroform pointing towards 
the center of the aromatic benzene ring. In order to con- 
nect to that study we also carry out calculations with the 
H atom pointing towards graphene and find that configu- 
ration to have a smaller (less favorable) adsorption energy 
than configurations with the H atom pointing away from 
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graphene. In this work we therefore focus on adsorption 
configurations with CI (or Br) atoms closest to graphene 
as shown in Figure [T] 

We define the adsorption energy E a as the difference in 
total energy when the molecule is adsorbed on graphene 
(-^ads fix-cell) an d when it is in the gas phase far away 
from graphene (E^ fix _ ccU ) 

_ P 77>t0t TTltot (T\ 

— -C'ads, fix-cell ^gas, fix-cell- I 1 J 

Here we follow the sign convention that yields positive 
values of E a for systems that bind. The two total-energy 
terms in ([!]) are both calculated with the adsorbant and 
graphene within one unit cell, the unit cell having the 
same size in both calculations. 1 Because E^°l^ flx _ ccll is 
calculated with the same adsorbant-adsorbant separation 
as ^ads fix-ceil tne direct lateral interaction is subtracted 
from our results [THl [T7] . 

The optimal positions of the atoms within chloro- 
form (and bromoform) are determined by minimization 
of the Hellmann-Feynman forces acting on the chlo- 
roform (or bromoform) atoms, when adsorbed on to 
graphene ("ads, fix-ceil") and when away from graphene 
( "gas, fix-ceil" ) . We use the molecular-dynamics optimiza- 
tion method "fast inertial relaxation engine" (FIRE) [TH] 
with the requirement that the size of the remaining force 
on each atom is less than 0.01 eV/A. The positions of 
the graphene atoms are left unchanged. The Hellmann- 
Feynman forces are derived from gradients in the electron 
density n. The optimization yields the bond lengths and 
angles in the molecules (after adsorption and in the gas 
phase) and the optimal position of the molecule with re- 
spect to graphene. The potential well for the molecule 
near graphene is very shallow. We therefore start the op- 
timization process in several different lateral and vertical 
positions of chloroform. 

We use an orthorhombic unit cell of size 3^3 a g x 3 a g x 
23.0 A for most of our calculations, but we also test the 
adsorption energy convergence with unit cell size, using a 
unit cell of size 3y/3a g x 5 a g x 23.0 A. Here a g — V% a U; 
and ao = 1.43 A is the clean graphene lattice constant as 
found earlier by relaxing the lateral size of the unit cell 
|13j . In the direction perpendicular to graphene there 
is ~ 12 A of vacuum above chloroform when desorbed 
(and ~ 19 A when adsorbed) to avoid vertical interaction 
across unit cell boundaries. 

The (valence) electron density and wave functions are 
represented on evenly distributed grids in space. To en- 




FIG. 1: Illustration of chloroform adsorbed on graphene for 
the 3^ a g x 3 a g unit cell (top panel) and the 3\/3 a g x 5 a 9 
unit cell (middle and bottom panels). The unit cell is outlined 
by the broken lines. Also shown are atoms of part of the 
neighboring repeated unit cells, illustrating the separation of 
the repeated images of the chloroform molecule. C atoms are 
represented by medium size gray circles, H atoms by small 
white circles, and CI atoms by large green circles. 



In previous work the change of the adsorbant structure from the 
deformed structure (after removal from graphene) into the gas 
phase structure was sometimes 13 calculated with a GGA ap- 
proximation of E xc because that approximation is less sensitive 
to changes in grid positions |13H15| . However, with the fine real- 
space grid used here and the rather small adsorbant molecules 
the total energy difference between the deformed and the gas 
phase structure of chloroform or bromoform is less than 1 meV. 
We therefore here solely use vdW-DF in the calculations. 



sure a good accuracy in our calculations we choose the 
density grid to have approximately 0.12 A grid point sep- 
aration in all three directions [19] . In all calculations we 
use a 2x2x1 Monkhorst-Pack A-point sampling of the 
Brillouin zone. 

In DFT the total energy E tot [n] is given as a functional 
of the electron density n. DFT is in principle exact, but 
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in practise the exchange-correlation part E xc of the total 
energy must be approximated. The exchange-correlation 
energy may also be written as a sum of exchange and 
correlation contributions, E xc = E x + E c . 

In the vdW-DF method the E x is chosen as the ex- 
change part of a generalized gradient approximation 
(GGA). In the original version of vdW-DF 5J (which 
we use here) the exchange chosen is that of the revPBE 
approximation [20] . The correlation energy E c is split 
[2"T] into a nearly-local part E® and a part that includes 
the most nonlocal interactions E^r, 

E c [n]=E° c [n}+Ef[n}. (2) 

In a homogeneous system the term E® is the correla- 
tion DA obtained from the local density approximation 
(LDA), and in general [5] we approximate E® by £' l i jDA . 
The term 

Ef[n] = ^J J drdr' n(r)0(r, r')n(r') (3) 

describes the dispersion interaction and vanishes for a 
homogeneous system. It is given by a kernel <fi which is 
explicitly stated in Ref. 5, 

The term E^ 1 is sensitive to changes in the local real 
space grid, for example when the nuclei positions are 
translated by a distance not corresponding to an integer 
number of real-space grid points [TH HHJ For cal- 

culations involving small (few meV) energy differences, 
like the calculations for the potential energy surface in- 
troduced in Section |IIID| we therefore keep positional 
changes of the rigid molecules to an integer number of 
grid points. 

Finally, we make sure that each of the total-energy 
calculations are accurately converged to a change in the 
total energy of less than 0.1 meV per unit cell, or less 
than approximately 10 -6 eV per atom in the unit cell, 
in the last three DFT iterations. This choice of allowed 
change is significantly smaller than the default of GPAW 
and is essential in order to discuss energy changes in this 
adsorption system where relevant differences in total en- 
ergies are down to the meV scale. 

III. RESULTS AND DISCUSSION 

Besides the total energies, in our calculations we also 
determine the molecular structure of the chloroform and 
bromoform molecules, both in the dilute gas phase and 
in the adsorbed phase. We compare the bond lengths 
and angles of the gas phase molecules with experimental 
values, and we determine the changes that occur when 
the molecules are adsorbed. 

The chloroform and bromoform molecules may adsorb 
in various orientations and positions on graphene. We de- 
termine the optimal orientation and adsorption distance, 
and discuss the influence of the adsorption position on 
the adsorption energy. Further, we test the convergence 



of the adsorption energy ([I]) with respect to lateral size 
of unit cell, i.e., the length of the smallest molecule-to- 
molecule distance. 



A. Structure of the desorbed and adsorbed 
molecules 

For the desorbed molecules the atomic positions within 
the molecule are determined in the process of determin- 
ing the total energy of the system of a molecule far away 
from graphene. The total energy of these gas phase struc- 
tures, along with the clean graphene total energy, are 
used as the desorbed system energy with which the ad- 
sorbed system total energies are compared. The sum of 
total energies of the desorbed system is used as the zero 
point of the adsorption energy. In practise for our des- 
orption calculations we keep graphene and chloroform (or 
bromoform) within the same unit cell, but far apart. 

From the atomic positions in the gas phase the bond 
lengths and bond angles within the chloroform and bro- 
moform molecules are extracted and listed in Table [J 
We find that the C-Cl bonds are slightly shorter than 
the C-Br bonds, a result that is expected because Br is 
a larger halogen atom than CI. The bond values we find 
are in reasonable agreement with the experimental val- 
ues listed in the NIST database 23]. All bond lengths 
are within 1-3% of the experimental values. Like for the 
experiments, the bond angles ZC1-C-C1 (and ZBr-C-Br) 
are a few degrees larger than the ZH-C-C1 (and ZH-C- 
Br) angles. For the bond angles our results deviate less 
than 1% from experiment (Table |l|. 

When adsorbed, the relative positions of the atoms in 
the adsorbant change. However, as we have also found in 
other small physisorbed molecules [T21 [H] , the changes 
are very small. Table [I] lists the changes. 

B. Adsorption energies 

Our main calculations are optimizations of the ad- 
sorbed chloroform molecule when it is positioned with 
the H atom pointing away from graphene ("H up"), as 
illustrated in Figure[T] We find adsorption energies in the 
350-360 meV range depending on the precise position on 
graphene. For the position in Figure [l] we find E a = 357 
meV. As discussed below, the differences in adsorption 
energies in the various positions on graphene are small. 

Like the structural changes in the chloroform and bro- 
moform molecule upon adsorption, the energetic changes 
are small. In fact, we find the contribution to the adsorp- 
tion energy from the deformation of the adsorbant to be 
about 0.2 meV, less than the accuracy of our results (~ 1 
meV). 

Table [TTJ also lists the adsorption energies obtained 
through experimental measurements. These find that 
chloroform binds stronger to graphene than our results, 
with a deviation of our results from the experiments 36% 
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TABLE I: Bond lengths and bond angles of the gas phase and the adsorbed phase chloroform and bromoform molecules. 
Adsorption structures are for the "H-up" structure (see text). For the theory results the vdW-DF method is used [5j [6]. Error 
bars on results from the NIST data base [23! are not available for chloroform. 



NIST This work 
gas phase gas phase ads. phase |A(gas — ads.) A(NIST — gas)| 


chloroform 












d C -H (A) 


1.073 


1.087 


1.087 


< 0.001 


0.01 (1%) 


dc-ci (A) 


1.762 


1.798 


1.801 


0.003 


0.03 (2%) 


ZC1-C-C1 (deg) 


110.92 


111.4 


111.1 


0.3 


0.5 (< 1%) 


ZH-C-C1 (deg) 


107.98 


107.7 


107.7 


< 0.1 


0.3 (< 1%) 


bromoform 












d C -H (A) 


1.11 ± 0.05 


1.085 


1.085 


< 0.001 


0.02 (2%) 


dC-Br (A) 


1.924 ± 0.005 


1.974 


1.977 


0.003 


0.05 (3%) 


ZBr-C-Br (deg) 


111.7±0.4 


111.5 


111.8 


0.3 


0.2 (< 1%) 


ZH-C-Br (deg) 


107.2 ± 0.4 


106.9 


106.9 


< 0.1 


0.3 (< 1%) 



TABLE II: Adsorption energies E a from theory and experiment, distance of molecular C atom from the plane of graphene, 
d,R, adsorption configuration (H atom sticking up or down), and unit cell used in calculations. The method vdW-DF of Refs. 
[U[6] is used when not noted otherwise. The adsorption structures for the experiment results are not known to us and for each 
molecule experiments are entered in the first entry line. We use orthogonal unit cells and a graphite lattice vector a g = V3ao 
with ao = 1.43 A. 

This work Experiments 



Structure Unitcell E a cLr E a 











[kJ/mol] 


[meV] 


[A] 


[kJ/mol] 


chloroform H up 


3^3 


X 


3 


34.4 


357 


4.20 


54±3 a , 36 A b 


H up 


3^3 


X 


5 


34.3 


356 


4.20 




H down 


3\/3 


X 


3 


33.5 


347 


3.48 




bromoform H up 


3^3 


X 


3 


39.0 


404 


4.32 




methane c 


3^3 


X 


3 


14.6 


152 


3.64 


13. 6 d , 17±l a 



"Thermal desorption spectroscopy measurements (at one monolayer), Ref. [8|. 
''Single atom gas chromotography, Ref. [9]. 

c Theory results from previous vdW-DF study with the same settings, Ref. [13] . 
d Temperature programed desorption (results extrapolated to isolated adsorbant), Ref. [25] . 



for thermal desorption spectroscopy [5] or 5% for single 
atom chromotography [9]. 

For our calculations we also report in Table [IT] the dis- 
tance dpi between the C atom of chloroform (or bromo- 
form) and graphene, in the adsorbed configuration. 

The majority of our calculations are carried out in a 
unit cell of lateral size a g x3a g , leading to a smallest 
(center-of-mass) molecule-to-molecule distance 7.43 A for 
the periodically repeated images, or the distance 5.00 A 
between two closest CI atoms each residing on a differ- 
ent chloroform molecule. In our calculational procedure 
described by (TJ) we subtract any direct interactions be- 
tween the molecules because both terms are calculated 
with the same lateral distance to neighboring molecules 
(i.e., with the same lateral size of the unit cell). Indi- 
rect interactions leading to changes in adsorption energy, 



could for example arise via a possible small deformation 
of the electron distribution on graphene, or via a possible 
but tiny deformation of atomic structure on chloroform. 
In order to check such possible effect we also calculated 
the adsorption of chloroform in a unit cell with 12.38 A 
molecule-to-molecule distance, the 3-\/3 a g x5a g unit cell. 
As evident from the results in Table \U\ there is very little 
difference between results of the small and large unit cell, 
confirming that the 3\/3a g x 3a g unit cell size is suffi- 
cient for size-convergence, provided that direct molecule- 
molecule interactions across unit cell boundaries are can- 
celled like we do here. 

In order to check whether the chloroform configura- 
tions with the H atom pointing to graphene are ener- 
getically more favorable than the "H up" configurations 
we also carry out a number of "H down" calculations. 
We calculate E a for the configuration with the H atom 
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above the center of an aromatic ring in graphene and in a 
number of configurations with the H atom near or above 
a graphene C-C bridge. We find the configuration with 
H above the center of an aromatic ring to be the most 
favorable of the H-down configurations (H centered con- 
figuration listed in Table [n| by up to 12 meV, but less 
favorable that configurations with the H atom pointing 
away from graphene. 

Previously, the dimer of chloroform with benzene has 
been studied using vdW-DF [7J. In that study the in- 
teraction between an aromatic 7r-system (represented by 
benzene) and an aliphatic C-H group (in chloroform) was 
at focus, and accordingly the orientation of the chloro- 
form was chosen to have the H-atom pointing towards the 
(center of) the benzene molecule. In Rcf. [7\ the binding 
energy 5.11 kcal/mol (21.4 kJ/mol or 222 meV/dimer) 
was found, with the distance from chloroform-C to the 
benzene-plane da = 3.6 A. The vdW-DF results of that 
study were compared with coupled cluster [CCSD(T)] 
calculations [26] with binding energy 5.60 kcal/mol (23.4 
kJ/mol or 243 meV/dimer) at dn = 3.2 A. In our cal- 
culations the interaction of chloroform is with the full 
graphene plane (as far as the vdW forces reach) and not 
only with the benzene molecule, we should therefore ex- 
pect a larger interaction energy than for the molecular 
dimer. Indeed, with an adsorption energy of 347 meV 
(Table [TTJ) we do find stronger binding than in the dimer 
case, stronger by about 123 meV, and with an adsorption 
position closer to the aromatic (graphene or benzene) 
ring by about 0.1 A compared to the vdW-DF calcula- 
tion and further away by approximately 0.3 A compared 
to the CCSD(T) calculation. 

In Ref. [4] the adsorption of chloroform on to the 
(5,0) and (8,8) single walled carbon nanotubes (SWCNT) 
was addressed with the use of the LDA approximation 
to DFT. Adsorption energies 200 meV [for the (5,0) 
SWCNT] and 150 meV [for the (8,8) SWCNT] were 
found. However, even though LDA seems to bind vdW 
materials it cannot be used for the inclusion of vdW in- 
teractions, as already pointed out by Harris |27) and dis- 
cussed also in Refs. [TBI [28]. The LDA results of Ref. @] 
are therefore not further discussed here. 



C. Adsorption of bromoform 

The Br atom is similar to CI but is heavier, i.e., has 
more electrons and larger molecular polarizability. We 
therefore expect a stronger binding to graphene than for 
chloroform. Indeed, our results in Table [TTJ show the ad- 
sorption energy to be 47 meV larger than that for chlo- 
roform. The adsorption distance as measured from the 
bromoform C atom is larger, reflecting the fact that the 
Br atom has a larger volume. The same trend, in the op- 
posite direction, is seen for methane [13] : the H atom has 
less electrons than CI, the binding is less strong (by ^200 
meV compared to chloroform) and the binding distance 
da is smaller because H has less volume than CI. 
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FIG. 2: Potential energy surface (PES) for chloroform on 
graphene in the "H-up" adsorption structure. Chloroform is 
here kept at a distance d — 4.20 A from graphene, measured 
from the chloroform C atom to the plane of graphene. The 
PES plot scans 1/4 of the 3v3a 9 x 3a g unit cell, with origo 
taken as the chloroform position of the top panel in Figure 
[TJ and with the same orientation as that panel. The energy 
scale measures the deviation from the global minimum. 



For the bromoform result shown in Table [TTJ we used 
the same position on graphene as illustrated in the top 
panel of Figure [TJ as the starting point for the structural 
relaxations. 



D. Chloroform potential energy surface 

In Figure[2]we show the potential energy surface (PES) 
for chloroform. The figure is obtained from translating 
chloroform along graphene with the chloroform C atom 
kept at a fixed distance d — 4.20 A from the graphene 
plane. In Figure [2] origo corresponds to the position 
shown in the top panel of Figure [TJ We find that the vari- 
ation in adsorption energy along graphene is small, with 
an energy difference of less than 6 meV when d = 4.20 A. 
This illustrates that the corrugation of graphene, as ex- 
perienced by the adsorbed chloroform molecules, is very 
small, and it takes only very little kinetic energy to over- 
come the barriers for lateral motion on graphene. There- 
fore, the concept of "adsorption sites" is not relevant in 
these physisorption studies [29, 30 . In effect, the chloro- 
form is free to move along graphene at all temperatures 
relevant in practical applications mentioned in the intro- 
duction. 

From Figure [2] we also find that the position used for 
the calculated adsorption energy in Table [TTJ is close to 
but not quite the optimal position, albeit not the ener- 
getically worst position either. In any case, the effect of 
our choice of position is a few meV, or less than 1% of 
the adsorption energy. 

In all calculations for the PES the relative positions 
of the atoms within the adsorbants are kept the same, 
thus no further relaxation of the molecular structure is 
allowed. However, since the change in molecular struc- 
ture when desorbing is on the sub-meV level, and less 
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than the accuracy of our calculations, we expect the dif- 
ferences in molecular structure in the various adsorption 
positions to be even smaller and the effect on the adsorp- 
tion energy to thus be negligible. 

All data points for the PES are from adsorbant posi- 
tions translated an integer number of grid points along 
the surface, that is, on a uniform orthogonal grid with 
0.12 A between grid lines. We restrict the possible posi- 
tions of chloroform in order to keep the atomic positions 
relative to the density grid the same in all calculations, 
thus avoiding any effects of the sensitivity of the vdW-DF 
results on the positioning on the grid. This is potentially 
important because the PES is mapping a very shallow 
energy landscape. 

The data for the PES are calculated for all grid po- 
sitions within a y/3 a g x 1 a g part of the 3\/3 a g x 3 a g 
unit cell, that is, by scanning over twice the area of the 
graphene primitive cell (which has only two C atoms). 
Therefore, when translating the chloroform molecule over 
the \/3 a g x 1 a g area we calculate two data points for each 
unique adsorption position. To lower the sub-meV noise 
we use the average of the total energies in the two equiv- 
alent positions for the plot in Figure [2] For clarity we 
also include in Figure [2] repetitions of the calculations in 
the two lateral directions. 

In Figure [2] only results of the fixed distance d = 4.20 
A are shown, but our calculations include also distances 
4.20±0.12 A in the same lateral positions as used in Fig- 
ure [2] This corresponds to moving the molecule one full 
grid spacing closer to or further away from graphene and 
redoing the PAS calculations. However, all adsorption 
energies at those distances are smaller than any of the 
adsorption energies at the d = 4.20 A distance, and it is 
clear that the optimal distance dn from graphene varies 
much less than 0.12 A when moving along graphene. Al- 
though the vertical part of the adsorption potential is 
shallow compared to covalent and ionic binding the lat- 
eral part is even more shallow. 

The differences in smallest and largest total energy 
within each of the three chloroform- C-to-graphene dis- 
tances are 7.7 meV, 4.8 meV, and 2.9 meV for the dis- 
tances d = 4.08 A, 4.20 A, and 4.32 A, respectively. Thus 
the corrugation of graphene, as seen by chloroform, be- 
comes slightly more pronounced in positions closer to 
graphene, even though the corrugation is small at all 
three distances considered here. 



E. Implications for environmental research. 

Even water from untreated drinking water wells con- 
tains THMs [3]. The THMs are spread in the environ- 
ment since chlorinated water is used for watering, and it 
leaks from swimming pools and enters waste water. It 
is also produced by salt used on winter roads. Chloro- 
form is also relatively volatile and escapes from water 
into the air with vapor (where inhalation may pose a 
health problem) or enters through the skin, for exam- 



ple during showers or in swimming pools. Absorption 
through the gastrointestinal tracts is fast and extensive, 
with the majority of ingested chloroform recovered in ex- 
pired air within a few hours [2 . With chloroform found 
in 11.4% of public wells in the U.S. [3J and chlorination 
still an important candidate for improving the quality of 
drinking water in developing countries, the occurrence of 
THMs is a potential human health concern, and methods 
to remove THMs after chlorination and before use of the 
water should be improved. 

Chlorination of water gives rise to a number of THMs 
as byproducts, mainly CI- and Br-based THMs. In or- 
der for these to be removed from water by adsorption 
on to graphene it is necessary that the adsorption en- 
ergy at least exceeds that of water on graphene, and that 
the adsorption energy is much higher than the barrier for 
thermal desorption at relevant temperatures, here around 
300 K which corresponds to -26 meV. In a CCSD(T) 
study [3T] the adsorption energy of water on graphene 
was found to be 135 meV, and in a recent vdW-DF study 
[3 2) (utilizing a different exchange functional than here, 
the optB86b exchange [33]) the water adsorption energy 
was found to be 140 meV. Both results are clearly smaller 
than our chloroform and bromoform adsorption energies 
that are in the range 350-410 meV. Thus, while not di- 
rectly giving proof that graphene can be used for water 
filters after chlorination our results do suggest that this 
may be possible. 



IV. SUMMARY 

We present a study of the adsorption on graphene of 
the two THMs chloroform and bromoform using the van 
der Waals density functional method vdW-DF. We find 
that chloroform and bromoform physisorb with their H 
atom pointing away from graphene, yielding adsorption 
energies 357 meV (34.4 kJ/mol) for chloroform and 404 
meV (39.0 kJ/mol) for bromoform. This suggests that 
these THMs bind sufficiently strongly to graphene for 
graphene to be used in filtering of chlorinated water to 
remove the THM byproducts. 
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